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An evidence that the dielectric aging in the polydomain Pb(Zr .45Ti . 55)03 thin films is controlled by pro- 
gressive pinning of 180° domain walls is presented. To provide such a conclusion, we use a general method, 
which is based on the study of the time evolution of the nonlinear, but anhystcrctic, dielectric response of the 
ferroelectric to a weak electric field. A thermodynamic model of the ferroelectric system where the dielectric 
response is controlled by bending movements of pinned 180° domain walls is developed. Within this model, 
the nonlinear permittivity of the ferroelectric is expressed as a function of the microstructural parameters of 
the domain pattern. It is shown that by using the analysis of the time evolution of the nonlinear permittivity, 
it is possible to estimate changes in the concentration of the pinning centers that block the movements of the 
180° domain walls during aging in polydomain pcrovskite ferroelectrics. 
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I. INTRODUCTION 

Aging of dielectric properties represents an unwanted 
feature in ferroelectric ceramics and lead zirconatc- 
titanatc (PZT) thin films, in particular, since it prevents 
their use in devices, which require very stable functional 
properties. It is believed that the origin of the dielectric 
aging in PZT and other similar perovskite ferroelectrics is 
caused by a rearrangement of pinning centers that block 
the movements of domain walls at more points and re- 
duces the domain wall (extrinsic) contribution to per- 
mittivity of the ferroelectric. 1 In the following text, we 
will call this concept aging due to progressive pinning. 
Unfortunately, examples of credible evidence for the pro- 
gressive pinning concept are still rather scarce 1 ' 2 due to 
strongly limited possibilities of experimental techniques. 
Even the direct observation of domain wall pinning itself 
is very rare and requires quite advanced experimental 
techniques 3 ' 4 . For this reason, progress in understanding 
the detailed role of domain wall pinning in phenomena 
such as fatigue, aging, or imprint 5 is rather complicated, 
since it relies mainly on indirect measurements. One pos- 
sibility to mention here is the study of the domain wall 
dynamics using measurements of the small-signal dielec- 
tric response, 6-8 which is very sensitive to the pinning of 
the domain walls. Unfortunately, these experiments were 
analyzed using models, which cannot serve any quantita- 



tive information on concentration of the pinning centers 
at the domain wall, etc. 

The aforementioned issues have motivated the study 
presented below, where we will apply a general method 
introduced recently in Ref. 9 that allows us to provide 
evidence that the dielectric aging of polydomain PZT 
films is controlled by progressive pinning of 180° domain 
walls. The adopted method is based on the measure- 
ment of the nonlinear dielectric response of non-polar 
polydomain ferroelectric samples to a weak electric field. 
This approach has several advantages. First, the appli- 
cation of a weak electric field to the ferroelectric sample 
has a minimum effect on the aging process. It is known 
that the electric cycling of the ferroelectric sample with 
electric fields comparable (in magnitude) to the coercive 
field has a strong dcaging effect. 10 Second, the dielec- 
tric response of the polydomain ferroelectric sample to 
the weak electric field is anhysteretic, 11 since it is con- 
trolled only by a fast reversible movement of the domain 
walls. This makes it possible to adopt a much simpler 
theoretical treatment in our analysis. It follows from the 
symmetry reasons that in the limit of the electric field 
tending zero, the permittivity ej of the polydomain non- 
polar ferroelectric sample (i.e., with the equal volumes 
of domains with the vectors of spontaneous polarization 
oriented along and against the applied electric field) is 
quadratically dependent on the electric field E. 
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where El is the small-signal permittivity and b is the di- 
electric nonlinearity constant. This is in contrast to the 
essentially hysteretic dielectric response of the ferroelec- 
tric to the subswitching electric field, which is accom- 
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panied by the Raylcigh-type linear electric field depen- 
dence of permittivity due to the irreversible movement 
of the domain walls. 12-17 Nevertheless, one should note 
that there may arise experimental situations where even 
the practically small applied electric field breaks the lim- 
its for the irreversible movements of the domain walls 
and, therefore, violates the condition of the limit of the 
electric field tending zero. Thus, one should very care- 
fully check that the experimental conditions satisfy the 
conditions for the applicability of the model. 

The final advantage of the adopted method to study 
the nonlinear permittivity of the ferroelectric polydomain 
system in a weak electric field is that in this case, we 
do not need to strictly identify the microscopic origin of 
the domain wall pinning. The reason is that in weak 
electric fields the extrinsic contribution to the permittiv- 
ity is controlled by fast reversible bending movements of 
the domain walls disregarding the nature of the pinning 
mechanism. Even if, in a particular sample, the domain 
wall pinning is strong due to many isolated crystal lattice 
impurities or it is weak due to the fluctuations of random 
fields in the vicinity of the domain wall, in a weak elec- 
tric field, the pinning effect of both the random bonds or 
random fields results only in the bending of the domain 
wall. It means that our analysis can be reduced down 
to the problem of the identification of the bending move- 
ments of pinned domain walls. For that reason, the key 
elements of our analysis is the development of the model 
for the description of bending movements pinned 180° 
domain walls. Within this model, we consider that the 
progressive pinning changes the bending condition of the 
domain wall during aging, which affects both the small 
signal permittivity El and the dielectric nonlincarity con- 
stant b. 

We will show that if the dielectric aging is caused by 
progressive pinning and is controlled by bending move- 
ments of 180° domain walls, there exists the following re- 
lation between time dependencies of the small-signal per- 
mittivity £h(t) and the dielectric nonlinearity constant 
b(t): 

y/bffi<xe L (t)-e Ci (2) 

where e c is the time-independent permittivity of the crys- 
tal lattice along the ferroelectric axis. To achieve the re- 
sult announced above and given by Eqs. (1) and (2), we 
first present, in Sec. II, the details of our model of a ferro- 
electric film, where the movements of 180° domain walls 
are locally blocked by pinning centers. By applying a 
straightforward thermodynamic methodology presented 
in Sec. Ill, we calculate the linear and nonlinear parts of 
the extrinsic contribution to permittivity controlled by 
bending movements of the pinned 180° domain walls. In 
Sec. IV, we show that in the progressive pinning aging 
scenario there exists a characteristic relation between the 
linear and nonlinear domain wall contributions to permit- 
tivity of the polydomain ferroelectric. We demonstrate 
that by using measurements of the time evolution of the 
nonlinear permittivity it is possible to obtain evidence 




FIG. 1. Model of 180° domain wall bending. Planar po- 
sition of the domain wall (doted line) is locally blocked by 
pinning centers (gray circles). Orientation of the vector of 
spontaneous polarization is indicated by gray arrows. When 
the external electric field E is applied to the ferroelectric film, 
the domain wall is bent between pinning centers and the pro- 
file of the wall deflection is described by the function u(x, y). 



on whether the dielectric aging is controlled by the pro- 
gressive pinning of 180° domain walls or not. Section V 
presents an application of our theory to aging experi- 
ments in [lll]-oriented PZT (45/55) thin films. 



II. MODEL OF 180° DOMAIN WALL BENDING 

Figure 1 shows the model of a ferroelectric film with 
a lamellar 180° ferroelectric domain pattern, where the 
average distance between the domain walls (domain wall 
spacing) is denoted by the symbol a. We consider that 
the top and bottom electrodes of the film are perpendic- 
ular to the ferroelectric axis x of the attached Cartesian 
coordinate system. Using the "hard ferroelectric" ap- 
proximation, we express the electric displacement within 
each ferroelectric domain as a sum of the spontaneous po- 
larization Pq (whose orientation differs from domain to 
domain) and the linear dielectric response of the crystal 
lattice to the electric field, 

D x = e e c E x ±P , (3a) 
D y>z — sos a Ey yZ , (3b) 

where e c and e a are the components of the permittivity 
tensor of the crystal lattice in the directions along and 
perpendicularly to the ferroelectric axis, respectively, and 
Sq is the permittivity of vacuum. 

In the absence of the electric field, we consider the do- 
main pattern to be neutral, i.e., that the volume fractions 
u+ and V- of the adjacent antiparallel domains are the 
same, and that the 180° domain walls are parallel to the 
x-y plane and pinned by immobile pinning centers, which 
are distributed in the x and y directions with an average 
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distance r between them. Thus, the pinning center den- 
sity on the domain wall is proportional to 1/ r 2 . 

When the electric field E is applied to the ferroelec- 
tric film, it exerts a pressure —2PqE on the domain wall 
and it bends the originally planar domain wall between 
the pinning centers. The domain wall bending produces 
a change in the volume fractions v + and t>_ of the do- 
mains with the spontaneous polarization oriented along 
and against the applied electric field, respectively. On 
the other hand, the domain wall deflection produces an 
increase in the domain wall area. Since the bent domain 
wall is no longer parallel to the vector of spontaneous 
polarization, a bound charge <7b = — APo,„ appears on 
the domain wall due to the discontinuous change in the 
normal component of spontaneous polarization at the do- 
main wall. If we denote the domain wall displacement in 
the z direction by the function u(x, y), its exact profile 
is given by minimizing the thermodynamic potential G w , 
which consists of the depolarizing field energy, the en- 
ergy associated with the crystal lattice polarization, the 
energy of the domain wall, and the energy supplied to 
the system by the external electric source. In this model, 
it yields the solution of Euler-Lagrange equation, which 
was analyzed earlier. 3 

In order to calculate the extrinsic contribution to the 
permittivity due to bending movements of pinned domain 
walls, we approximate the generally random distribution 
of pinning centers by a periodic one with equidistant pin- 
ning centers, but with the same average density 1/r 2 . We 
consider that the pinning centers divide the infinite area 
of the domain wall into square segments with edges par- 
allel to the x and y directions in such a way that each 
segment of the domain wall is pinned at its corners and 
the domain wall can move freely in the interior of each 
segment, as indicated in Fig. 1. This approximation re- 
flects the assumption that the main contribution to the 
permittivity is coming from the structures on the do- 
main wall with the periods that are close to the mean 
period in the system r. The strict procedure would be 
to expand the random distribution of the pinning centers 
in the Fourier series. Our approximation is good if the 
Fourier spectrum is not very wide. The consideration 
of the periodicity of pinning centers in the direction of 
the film thickness also limits the applicability of the do- 
main wall bending model to samples of thicknesses much 
greater than the average pinning centers distance, i.e., 
h > r. 

Finally, in most ferroelectric materials as well as in the 
samples used in our experiments, it is easy to show that 
even small domain wall displacements (about a lattice 
constant) are sufficient to produce the observed values of 
extrinsic contributions to the permittivity. Therefore, it 
is fully justified to consider that the maximum deflection 
r\ of the domain wall is much smaller than the average 
distance between pinning centers r. In this case, if we 
take the origin of the coordinate system in the middle of 
each square segment, it is convenient to approximate the 
profile of the domain wall deflection in this segment by a 



parabolic function of the form, 



u(x, y) = r] 



1 



a x 



,2 1 



r 2 (l + a) 



(4) 



where r\ is the maximum deflection of the domain wall 
and the parameter a is introduced in order to take into 
account the anisotropy of the radius of curvature of the 
bent domain wall due to the strong depolarizing effect 
in the direction of the ferroelectric x axis. Coordinate 
values of x and y are running over the interval from —r/2 
to r/2. 



III. EXTRINSIC PERMITTIVITY 

When the alternating electric field is applied to the 
ferroelectric sample, the domain walls start to vibrate 
between the pinning centers, which results in a change 
in the volume fractions v + and u_ of the antiparallel do- 
mains, where v+ + u_ = 1. This represents the source 
of the extrinsic contribution to the permittivity of the 
ferroelectric film. To calculate the extrinsic contribution 
to the permittivity, it is convenient to consider the net 
spontaneous polarization Pjv , which is given by the differ- 
ence in the volume fraction of domains with the vector of 
spontaneous polarization oriented along and against the 
applied electric field, i.e., Pjv = («+ — v~) Pq, and which 
can be expressed in the form, 



Pn = 



2Po 

nr 2 



u(x,y) dA= —P . 

3a 



(5) 



where the integral is taken over the square segment 
A, while x and y are running from —r/2 to r/2. In 
the following text, we will use the above expression to 
measure the maximum deflection of the domain wall 
in terms of the net spontaneous polarization, i.e., r/ = 
(3a Pjv)/ (4 Po). The response of the net spontaneous po- 
larization with respect to the applied electric field E is 
controlled by the thermodynamic function G w per unit 
volume of the ferroelectric, 



G w = $ s (Pjv) + $dc P (Pv) - PnE, 
where the functions 



(6) 



* s = ^2 J \j l + u K x ' y) + u l( x , v) dA ( 7 ) 



and 



$do P = [ 7:07, (p d Jl + ul{x,y) + u 2 (x,y) dA (8) 

ar J A z v 

represent two contributions to the thermodynamic func- 
tion G w due to the increase in the domain wall area and 
due to the depolarizing field, respectively. In Eqs. (7) 
and (8), the integrals are taken over the square segment 
A where x and y are running from —r/2 to r/2, the sym- 
bol a w stands for the surface energy density associated 
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with the surface tension of the bent domain wall, and 
the symbol ifd stands for the electrostatic potential on 
the domain wall associated with the depolarizing field, 
which is produced by the bound charges ov 

The function <fr s can be calculated in a straightforward 
way by direct substitution of Eq. (4) into Eq. (7). The 
leading terms of $ s with respect to the net spontaneous 
polarization equal 



<I>.s 



3a<7,„ (1 



4r 2 P 2 (1 + a) 



P z - 



9a 3 a w (9 + 10a 2 + 9a 4 



80 P$ r 



4 _4 



(l + a) 



-P 



N- 



(9) 



G w is dominated by the lowest (quadratic) term. Un- 
der this consideration, we can find the minimum of the 
function G w with respect to a and the condition given 
by Eq. (13) yields 




(14) 



In most of the samples, it is reasonable to consider 
that the domain wall thickness a w is much smaller than 
the average distance between the pinning centers r, i.e., 
a w <C r, and, in this case, Eq. (14) can be further sim- 
plified as 



In order to express the function ^dep, it is necessary to 
calculate the spatial distribution of the electrostatic po- 
tential ifd in the vicinity of the bent domain wall. Since 
this goes beyond the main scope of this paper, the de- 
tailed calculations of the functions <fd and $dep are pre- 
sented in the Appendix A. Here we present only the 
result for the leading terms of the function $dcp with 
respect to the net spontaneous polarization, 



0.174aa 2 



dep 



e a r(l + a) 2 ^/e a e c 



N 



3a 2 (l + 1.43a 2 
2P 2 r 2 (l+a) 2 



-P 



N 



(10) 



Now, the thermodynamic potential G w per unit vol- 
ume of ferroelectric can be expressed in the form of a 
Taylor expansion with respect to P/v- Disregarding the 
constant term, this function reads as 



3aP^ 



(a 2 + l) 



4r 2 (l + a) : 
9a 3 P* 



P 2 



(9- 



7.2 ra 2 
10a 2 + 9a 4 ) 



80P 2 r 4 (l+a) 4 
ra 2 (2.33 + 3.32a 2 ) 



P 2 



PnE, 



(11) 



where the first term in each square bracket is due to the 
increase in the energy of the surface tension of the bent 
domain wall and the second term is due to the increase 
in the depolarizing field energy produced by the bound 
charges at the domain wall. 

In what follows, it will be convenient to express the 
surface energy density of the domain wall a w in the form, 



Pq 

6e e c ' 



(12) 



where the parameter a w is of the order of the domain 
wall thickness. The unknown value of the parameter a is 
determined from the condition 



1.4r 



(15) 



After substitution of Eqs. (12) and (15) into Eq. (11), we 
obtain the following form for the expansion of the func- 
tion G w with respect to the net spontaneous polarization: 



G w (Pn, E) 



8eoe c r 2 



P 2 - 



0.17 a 3 a w 
e e c r 4 P§ 



P^-PnE. (16) 



It should be noted that in the case of most of the high- 
quality samples, a m i n 1 and the thermodynamic po- 
tential G w is dominated by the surface tension, which 
can be readily seen after substituting Eqs. (12) and (15) 
into Eqs. (9) and (10). 

The response of the net spontaneous polarization P/v 
to the applied electric field E can be found from the 
condition for the minimum of the function G w , 



dG w {P N , E)/dP N = 



(17) 



and it can be expressed as a Taylor expansion with re- 
spect to the applied electric field, 



P N (E)* A -^E- 



172.8 r^elel £3 
a al, P? 



(18) 



Now we use the fact that the average electric displace- 
ment of the polydomain film along the ferroelectric axis 
Df is given by the sum of the linear dielectric response of 
the crystal lattice to the electric field and the net sponta- 
neous polarization, i.e., Df(E) = eoS c E + Pm{E). Later 
we will check that the dielectric nonlinearity of the whole 
system is dominated by the bending mechanism and, in 
this case, the field dependence of permittivity of the fer- 
roelectric polydomain film can be expressed in the form 



(19) 



where e c is the intrinsic permittivity along the sponta- 
neous polarization and 



dG w {P N ,E;a)/da = 0. 



(13) 



Later, it will be checked that in the case of a stable do- 
main pattern and for weak applied fields, the function 



e w {r, a) 
b(r, a) 



1\ 518.4 r A elel 
aalP 2 



(20) 
(21) 
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are the coefficients of the small-signal linear and the 
quadratic terms of the extrinsic contribution to permit- 
tivity with respect to the applied electric field, respec- 
tively. Here it should be noted that Eq. (18), is actually 
derived for the case of the dc field dependence of the net 
spontaneous polarization on the applied electric field. On 
the other hand, the experimental part of the paper ad- 
dresses the ac field amplitude dependence of the average 
(mean) permittivity. Therefore, the relations given by 
Eqs. (20)-(21) are expressed for the ac field amplitude 
dependence of the average (mean) permittivity and the 
coefficient b differs from the corresponding term in Eq. 
(18) by a factor 1/4. 



IV. PROGRESSIVE PINNING OF DOMAIN WALLS 
AND DIELECTRIC AGING 

Now it is seen that values of the small-signal extrinsic 
permittivity e w and of the dielectric nonlinearity coeffi- 
cient b are controlled by the domain pattern configura- 
tion, i.e., by the pinning center average distance r and 
by the domain spacing a. It is natural to expect that 
the parameters r and a can evolve with time, which may 
represent a source of aging of the dielectric response and, 
therefore, a source of the time evolution of the field de- 
pendence of permittivity Ef(E). However, one should 
first identify the thermodynamic force, which can drive 
the system to evolve with time and which can be respon- 
sible for the aging process. 

First let us focus on the possible role of the domain 
spacing a in the process of aging. In real systems, the 
domain pattern is usually controlled by the prehistory of 
the sample so that the domain spacing may be essentially 
different from its equilibrium value at the given temper- 
ature. It means that there clearly exists a thermody- 
namic force that drives the system to reach the state with 
the equilibrium domain spacing. This thermodynamic 
force actually originates from the competition between 
the energy of the domain walls, the electrostatic energy, 
and the surface energy at the interface between the fer- 
roelectric and the electrode. 18-21 Nevertheless, frequent 
observations of rather stable but essentially nonequilib- 
rium domain patterns 22,23 represent a clear indication 
that these energies are usually much weaker than the en- 
ergies involved in the pinning-depinning processes during 
the domain wall movement under the application of sub- 
switching electric fields. Therefore, it is not very reason- 
able to consider that the possible change in the domain 
spacing could be responsible for the aging of the dielectric 
response. 

The second scenario, which can be described within 
our domain wall bending model, is the progressive pin- 
ning of the domain walls. Although it has been already 
mentioned in Sec. I that one can apply the domain wall 
bending model on the weak-field nonlinear permittivity 
data disregarding the origin of pinning, one particular 
mechanism in perovskitc ferroclectrics to mention here 



is pinning by the orientation of dipolc defects. 1 This 
model is based on the interaction of the domain walls 
with the dipoles formed by an acceptor ion (e.g., Ni 2+ , 
Fe 2+ , etc.) at the Ti 4+ site and an oxygen vacancy in the 
surrounding oxygen octahedron. Since the free energy of 
the dipole defect in ferroelectrics depends on its orienta- 
tion with respect to the vector of spontaneous polariza- 
tion, the thermodynamic force — which drives the dipolc 
defects to align with the spontaneous polarization vector 
in each domain in order to minimize the free energy of 
the whole system — causes an increase in the number of 
pinning centers that blocks the domain wall movement 
under a weak electric field. Thus, the increase in the 
number of such pinning centers naturally results in the 
decrease in the average distance r between them. 

Our model for bending movements of 180° domain 
walls makes it possible to identify such aging process 
using measurements of nonlinear permittivity. We can 
express the average distance between pinning centers r 
from Eq. (20) and substitute it into the formula for the di- 
electric nonlinearity coefficient b given by Eq. (21). This 
gives the relation between the parameters b and e w in the 
form, 



p2 
1 II 



(22) 



If such a relation between e w and b, which is typical for 
this aging process and which is quite different from the 
prediction of, e.g., the Landau theory, is identified in the 
aging measurements, it gives a reasonable hint that the 
progressive pinning mechanism is responsible for the di- 
electric aging of polydomain ferroelectrics. Therefore, 
the principal result of this work is that by considering the 
domain wall bending mechanism and using the analysis 
of the time evolution of the nonlinear dielectric response, 
we can provide evidence that the dielectric aging is con- 
trolled by the progressive pinning of the 180° domain 
walls. 

In addition, we can study the evolution of the pinning 
centers concentration at the domain walls. By combin- 
ing Eqs. (20) and (21) and by eliminating the average 
domain spacing, one can estimate the time dependence 
of the average pinning center distance according to the 
following formula: 



-(<) 



0.175 a w P 



b{t) 



,{t) 



(23) 



V. EXPERIMENTAL VERIFICATION 

In this section, we present a direct way of the exper- 
imental verification of the considered model for bending 
movements of the pinned domain walls to PZT ferroelec- 
tric films of tetragonal symmetry. Our model can be 
directly applicable to [001]-oricntcd PZT films. Unfor- 
tunately, there exists a reason, which makes the use of 
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[001]-orientcd PZT films not very practical to demon- 
strate the relevance of the domain wall bending model. 
Namely, the domain composition of c and/or a domains 
is difficult to control in such films. Because of the large 
dielectric anisotropy of PZT, this makes it difficult to 
control the lattice contribution to the permittivity. On 
the other hand, it appears to be convenient to use [111]- 
oricnted films; since their lattice contribution is unique 
and independent on the domain pattern configuration, all 
180° domain walls have the same orientation with respect 
to the applied electric field and, thus, they contribute 
identically to the dielectric response. These features are 
important in applying our model to experimental data. 
Since the application of our model to [111] -oriented PZT 
films requires some modifications, we describe the ap- 
plication of the domain wall bending model to the both 
particular cases separately. 




1 2 
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FIG. 2. (Color online) ac-field dependence of the dielec- 
tric permittivity (real part) within four decades of time. The 
measuring frequency was 1 kHz. The scattered markers rep- 
resent the experimentally measured data. The solid lines are 
quadratic fittings using Eq. (19). 



A. Application to [001]-oriented PZT films 

To prove that the progressive aging scenario is respon- 
sible for the evolution of the dielectric response, we need 
to distinguish the extrinsic e w + bE 2 from the intrinsic e c 
contributions to the field dependence of the permittivity 
£f(E). Since it is natural to expect that the intrinsic con- 
tribution to the permittivity e c does not change in time, 
we can compare the measured small-signal permittivity 
£l, given by the formula, 



e L = £/(0) = e c + e % 



(24) 



with the dielectric nonlinearity constant b, which can be 
determined from the experimental data using the follow- 
ing expression: 



b =[e f (E) -e L ]/E 2 . 



(25) 



Then, with the use of Eq. (24), the relationship given in 
Eq. (22) can be rewritten in the form, 



8.1ege c / a \ 



(26) 



Finally, taking the square root of b from the above equa- 
tion, we arrive at the following relationship between el 
and b, which can be cross checked experimentally: 



Vb ~ Ks jj — B, 



where 



K 



■ l£o£c 

p2 

1 ii 



— , B = Ke c 



(27) 



(28) 



Therefore the validity of Eq. (22) can be demonstrated 
by a linear relationship between the values of Vo and the 
small-signal dielectric permittivity sl- 
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FIG. 3. (Color online) Linear dependence of the square root 
of the coefficient b* of the quadratic term in the expansion of 
the permittivity with respect to the applied electric field on 
the total weak-field permittivity of the ferroelectric film. 



B. Application to [lll]-oriented PZT films 

In the case of the [lll]-oriented PZT films, one can- 
not apply our model in a such a straightforward way as 
it was done in Sec. VA. There are two main reasons 
for that. First, the bending movements of pinned 180° 
domain walls contribute to the value of the permittivity 
in the direction of the vector of spontaneous polariza- 
tion. Second, perovskite ferroelectrics are materials with 
a rather large dielectric anisotropy in the directions par- 
allel e c and perpendicular e a to the orientation of the 
vector of spontaneous polarization. In order to obtain 
reasonable numerical estimation of the microstructural 
parameters of the domain pattern, one should take the 
aforementioned point into account while interpreting the 
dielectric nonlinearity measurements. In the following 
text, we will denote all physical quantities, which are 
measured on the [111] oriented film, by a star superscript 

The electric field dependence of the out-of-plane per- 
mittivity of the [111] -oriented PZT film in the tetragonal 
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phase e j (E) is given by the formula, 

s*(E) = E f (E) cos 2 {9) +e a sin 2 (6>) 



(29) 



= (e c + e w + bE 2 ) cos 2 {6) + e a sin 2 (0), 

where E is the electric field along the vector of spon- 
taneous polarization, the symbol 9 stands for the an- 
gle between the spontaneous polarization and the nor- 
mal to the plane of the film vector, which is equal to 
9 = arctan(i/2) = 54.7°. The value of E is given by 
the applied electric field in the [111] direction E*, i.e., 
E = E* cos 9. If we introduce the dielectric anisotropy 
factor £ by the formula 



£ = e a /e c , 



(30) 



the measured small-signal permittivity e* L and the dielec- 
tric nonlinearity constant b* are then equal to 



e* L = e w cos 2 (9) 
b* = b cos 4 (6>). 



cos 2 (6>) + £sin(0)' 



(31) 
(32) 



Now we can follow the same procedure as in Sec. VA, 
i.e., we express the value of e w from Eq. (31) and we 
combine it with Eqs. (22) and (32). Finally, we obtain 
again the linear relationship between \/b~* and e* L , 



K*e* L - B*, 



where 



K* 


= K, 


B* 


= e c K 


K 





cos(26>) 



1 (I 



(33) 

(34a) 
(34b) 

(34c) 



VI. RESULTS AND A DISCUSSION 

The presented model of 180° domain wall bending was 
applied to study the evolution of the dielectric response in 
a [lll]-oriented Pb(Zr .45Tio.55)0 3 thin film (240 nm in 
thickness) . The film was deposited via a standard sol-gel 
method on a Pt-coated Si substrate. 13 Desired platinum 
patterns were vaporized on the surface of the crystallized 
film to form the top electrode. The out-of-plane dielectric 
response was measured using a Hewlett-Packard (HP) 
4284A high-precision impedance analyzer. The film was 
first depoled with a fast-decayed low- frequency (1 Hz) ac 
field (the amplitude decays from 80 kV/cm to zero in five 
periods). The dielectric response was then recorded as a 
function of the ac driving field {Eq <3 kV/cm) and of 
the aging time at room temperature. The coercive field 
of the film is 60 kV/cm and, therefore, it is much larger 
than the maximum field used in this study. 

Figure 2 shows the ac-field dependence of the dielec- 
tric permittivity within four decades of time, which can 



be well fitted by the quadratic relation given in Eq. (19). 
This is in contrast to the Rayleigh-type relation, where 
the dielectric permittivity increases linearly with the ac 
field as a result of the irreversible movement of domain 
walls under subswitching fields (usually a few tens of 
kV/cm). 12-17 The bending of 180° domain walls is actu- 
ally a fast reversible process. Figure 3 shows the linear re- 
lationship between the values of \/b* and the small-signal 
dielectric permittivity e* L . Therefore, our experimental 
data unambiguously prove the validity of Eq. (33). 

Using the fit of results presented in Fig. 3, we can es- 
timate the lattice permittivity e c and the ratio a/a w . 
This is done in two steps. First, the slope K* = 
85 x lO^mV- 1 and the offset B* = 62 x lO^mV -1 of 
the linear dependence of Vb* versus e* h are determined 
using a standard linear regression analysis. Second, the 
system of Eqs. (34) is solved by taking the values of spon- 
taneous polarization Pq = 0.54 C/m 2 (Ref. 24 ) and the 
dielectric anisotropy factor £ = 3. 3. 25 Solution of the sys- 
tem of Eqs. (34) gives a value for the lattice permittivity 
e c equal to 280, which is in excellent agreement with the 
thermodynamic values, 24 ' 25 and for the ratio a/a w equal 
to 12 x 10 3 . From Eq. (31), the small-signal extrinsic con- 
tribution tothe permittivity along the spontaneous polar- 
ization e w is estimated to be decreasing from 388 down 
to 296 during the time of 10 4 s of the aging experiment. 
The bending mechanism contributes to the total weak- 
field permittivity along the spontaneous polarization by 
about 50%, which is also in agreement with the recent 
dielectric measurements of the weak-field permittivity of 
PZT ceramics in a very wide frequency range. 26 From 
the ratio a/a u , = 12 x 10 3 . we can estimate the average 
domain spacing using the domain wall thickness as found 
from the Landau-Ginzburg-Dcvonshire theory (0.7 nm) 
(Ref. 27 ) or by an ab-initio calculation (0.5 nm). 28 The av- 
erage distance between adjacent 180° domain walls thus 
has a value of about 5.9 fj,m. This specific distance is 
much larger than the typical grain size (~100 nm), indi- 
cating a low density of 180° domain walls. This is reason- 
able in the case of strongly clamped films, as was shown 
by the transmission electron microscopy (TEM) observa- 
tion in Ref. 14 In addition, the average distance between 
adjacent 180° domain walls is controlled by the prehis- 
tory of the sample and by the depoling process used be- 
fore the aging experiment. Finally, our experiments indi- 
cate that the average distance between the pinning cen- 
ters has decreased from 33 down to 29 nm during 10 4 s. 
These values, which are much smaller than the film thick- 
ness, correspond to the volume concentration of crystal 
lattice impurities being higher than 0.03%, which is ac- 
ceptable since it is known that the nominally pure PZT 
ceramic films possess naturally occurring acceptor impu- 
rities. 

Considering the key element of our method — the model 
of bending movements of pinned 180° domain walls — we 
have calculated the linear and nonlinear contributions to 
extrinsic permittivity in the polydomain ferroelectric and 
used this result to analyze our experimental data. As a 
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result, wc were able to extract some microstructural pa- 
rameters of the domain pattern, i.e., the average distance 
between the pinning centers and the domain spacing, and 
material parameters of the ferroelectric, i.e., the intrin- 
sic permittivity of the crystal lattice. However, three 
important assumptions have been made during the de- 
velopment of the theoretical model and now it should be 
checked whether the values of the parameters fitted from 
our experimental data do not violate such assumptions: 
(a) the nonlinearity is dominated by bending of the 180° 
domain walls, (b) the maximum deflection of the bent 
domain wall is much smaller than the average distance 
between the pinning centers, and (c) the free energy G w 
that controls the net spontaneous polarization response 
to the electric field Pn(E) is dominated by the lowest 
(quadratic) term. 

The first assumption can be verified by comparing 
the dielectric nonlinearity constant that is controlled by 
bending movements of the 180° domain walls b with the 
dielectric nonlinearity of the crystal lattice 6 C , which is 
equal to 9 



1 9 r 2 r 3 

° c — p2 



1 - 



3&» 
8e, 



By substituting the numerical results of our study, one 
can immediately see that b c = 1.6 x 10~ 16 m 2 V~ 2 , which 
is about 6 orders of magnitude smaller than the values 
observed in the experimental part of this study, i.e., b ss 
8 x l(r 10 m 2 V- 2 . 

The second assumption can bc verified by using 
Eqs. (5) and (18). After the substitution of the numerical 
parameters fitted from our experiments, we can obtain 
that, at the maximum fields applied to our sample, i.e., 
3 kV/cm, the maximum deflection rj of the domain wall is 
about 17 nm. This value is more than two times smaller 
than the estimated average distance between the pinning 
centers r (over 40 nm) . 

The third assumption can bc verified by the substitu- 
tion of Eq. (5) into Eq. (11). With use of the numeri- 
cal parameters fitted from our experiments, we can see 
that, at the maximum fields applied to our sample, i.e., 
3kV/cm, the value of the lowest (quadratic) term in the 
expansion of the thermodynamic function G w is about 
302 Jm" 3 . This value is more than three times greater 
than the value of the higher- (fourth-) order term, which 
is approximately -90 J m~ 3 . 

Therefore, it is seen that all assumptions made in the 
theoretical part of our study are fully justified and ap- 
plicable to real ferroelectric samples. By applying the 
main results of the theoretical model [Eqs. (19) to (22)] 
to our experimental data, we have provided a strong evi- 
dence that the dielectric nonlinearity in polydomain fcr- 
roelectrics is predominantly controlled by the considered 
bending mechanism. In addition, we have shown that 
from the evolution of the relation between the small- 
signal permittivity sl and the dielectric nonlinearity con- 
stant b it is possible to distinguish the microstructural 
mechanisms, which are responsible for the aging of the 



dielectric response. In our particular experimental case, 
the linear relationship between the values of Vb* and e* L 
provided us with evidence that the decrease in the linear 
dielectric permittivity during aging can be attributed to 
the increase in the average pinning center density on the 
domain wall, indicating a progressive pinning nature of 
the aging phenomenon. Finally, Eqs. (20) and (21) show 
that — from known values of the parameters el and b at 
a given time — it is possible to estimate the actual values 
of the domain spacing a and the average pinning center 
density on the domain wall 1/r 2 . Therefore, we believe 
that our results can be used as a simple and useful tool 
for getting a deeper insight into the configuration of the 
domain pattern and the quality of ferroelectric thin films 
and for providing a way to identify the evolution of the 
180° domain pattern microstructure in perovskitc ferro- 
electrics. 
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Appendix A: Electrostatic energy of the bent 180° domain 
wall 

In this appendix, we present the detailed calculation of 
the electrostatic potential <fd and the depolarizing field 
energy per unite volume of the ferroelectric $dcp with the 
bent 180° domain walls. When the domain wall is bent, 
the discontinuous change in the normal component of 
spontaneous polarization Po,n at the domain wall yields 
the appearance of a bound charge of surface density o\, = 
-AP 0l n, 



<7b{x, y) 



-2P u x (x 7 y) 



yjl + ul(x, y)+u 2 (x, y), 



(Al) 



where the functions u x (x,y) and u y (x,y) are the partial 
derivatives of the function u{x, y) with respect to x and y, 
respectively. If the maximum deflection rj of the domain 
wall is much smaller than the average distance between 
pinning centers r then by using Eq. (4), the bound charge 
density on the domain wall can be expressed as a Taylor 
expansion with respect to rj as follows: 



&b(x, y) 



\6P ax 512P ax(y 2 + a 2 x 2 ; 
r 2 (l + a) ?/ ~ r 6 (l + a) 3 



if. (A2) 



If we consider that the maximum deflection of the do- 
main wall is much smaller than the average distance be- 
tween the pinning centers, i.e., rj <C r, the electrostatic 
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potential ipd can be approximated by the solution of the 
electrostatic problem where we assume that the bound 
charges o\, on the bent domain wall are located at the 
original position of the domain wall, i.e., at z = 0. Then, 
the electrostatic potential <fd is given by the solution of 
the Laplace equation, 



d 2 ip d ( d 2 ip d d 2 (p d 



dx 2 



dy 2 



dz 2 



0. 



(A3a) 



with the internal boundary conditions for the continuity 
of the normal component of electric displacement and for 
the continuity of the electrostatic potential at the domain 
wall, i.e., at z = 0, 



d Vd +) _ d Vd ] _ (+) _ (-) 



dz 



£0£a 



(A3b) 



where the superscripts (+) and (— ) denote the electro- 
static potential for z > and z < 0, respectively. Consid- 
ering the periodicity of the bound charge surface density 
on the domain wall in both the x and y directions with 
the period equal to r, we will look for a solution of the 
system of Eq. (A3) in the form of a Fourier series. It can 
be easily shown that the functions 



¥><i (z, V, z) 



oo 

E < 

n,m—l 



) 



' 2n7r.T \ / 2miry 
1771 sin \ - I cos I — - — 



= =F(27rz/r) y / m 2 +n 2 e c /e a 



(A4) 



satisfy the Laplace (A3a). The unknown coefficients 4> n o 
and 4> nm can be found by substituting Eqs. (A2) and 
(A4) into Eq. (A3b). Straightforward calculations yield 
the formulae for the unknown coefficients, 



MO — 



3aa(-iy 



^■K 2 e 0y /£ a e c (l + a) 



(A5a) 



' 3a 2 [l8a 2 -n 2 7r 2 (l + 3a 2 )] n3 ) 
r'N H „ „ o „n — ijv > , 



2n 2 ir 2 r 2 P 2 (l + a) 2 
54a 3 a(-l) m+n 



m 2 n7r 4 e r 2 P ( 2 y / m 2 e 2 + n 2 e c e a (l + a) 3 



P 



N- 



(A5b) 

After substitution of Eq. (Al) into Eq. (8), the function 
<i>dcp can be expressed in the following forms: 

*do P = % [ u x (x,y)<pd(x, y, 0) dA (A6) 

ar z J A 

^ -6a(-l) n P N 



nnr (I + a) 

n—l x ' 



'MO- 



It should be noted that the coefficients 4> nm do not enter 
the formula for i»dep because the function u x (x,y) does 
not depend on y and the integration of all terms with 



cos(2rrnry / r) over its period r gives zero. After further 
substitution of Eq. (A5) into Eq. (A6) and performing 
the summation, the Taylor expansion of the $dep with 
respect to the net spontaneous polarization is equal to 



0.174aa 2 



dep 



e r(l + a) 2 ^/e a s c 



P z -■ 



3a 2 (l + 1.43a 2 
2P 2 r 2 (l + a) 2 



-P : 



N 



(A7) 
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